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THE THEORY OF PROPELLERS 
IV - THRUST, ENERGY, AND EFFICIENCY FORMULAS FOR SINGLE- 
AND DUAL-ROTATING PROPELLERS WITH IDEAL 
CIRCULATION DISTRIBUTION 
By Theodore Theodorsen 



SmiMARY 



Simple and exact expressions are given for the 
efficiency of single- and dual-rotating propellers with 
ideal circulation distribution as given by the Goldstein 
functions for single-rotating propellers and by the new 
functions for dual-rotating propellers from part I of the 
present series. The efficiency is shown to depend 
primarily on a defined load factor and, to a very small 
extent, on an axial loss factor. Tables and graphs are 
included for practical use of the results. The present 
r>aper is the fourth in a series on the theory of propellers 



INTRODUCTION 



The thrust, the energy loss, and the efficiency of 
a propeller are given completely and uniquely by the 
condition of the wake far behind the propeller. Detailed 
knowledge of the oropeller required to create the 
particular wake pattern is not needed; in fact, the pro- 
peller is not uniquely determined by the wake pattern. 
An element of lift may be transposed in a direction 
tangent to the vortex surface in such a manner as to 
maintain the identical vortex pattern far behind the 
propeller . 

Several equivalent propellers may thus exist - all 
with the same vortex surface far behind the propeller. 
Such quantities as the diameter, the pitch, and the rate 
of advance of the si;rface of discontinuity far behind 
the propeller are therefore of a more fundamental 
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significance in many respects than the similar quantities 
referring to the propeller. At any rate, it has been 
found convenient for the present theoretical treatment 
to consider all quantities as referring to the conditions 
of the ultimate wake. Only in the final stage of the 
actual design of the propeller are the interrelations 
of the propeller and the ultimate wake of concern. For 
the present investigation only knowledge of the ultimate 
wake is required. The thrust, the various energy losses, 
and the efficiency are dependent only on the ultimate 
wake . 

The present naper is the fourth in a series on the 
theory of propellers. The first of the series (refer- 
ence 1) deals with a set of new functions for the thrust 
distribution of dual-rotating propellers. The second 
(reference 2) concerns the axial interference velocity, 
and the third (reference 5) treats of the contraction 
of the propeller wake. 



SYMBOLS 



e 


angular coordinate on vortex sheet 


R 


tip radius of propeller 




X 


nondimensional radius in terms 


of tip radius 


z 


axial coordinate 




p 


number of blades of propeller; 


also, pressure 


r 


circulation at radius x 




00 


angular velocity of propeller 




V 


advance velocity of propeller 




w 


rearward displacement velocity 
surface (at infinity) 


of helical vortex 


\ 


. ^. /V 4- w\ 

advance ratio { i 




K(x) 


circulation function for single 
/ PFoo \ 
\2.Tr{ V + w )wy 
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K(x,0) circulation function for dual rotation 

K mass coefficient 

>1 . pi n2Tr 

K(x)x dx or i I K(x,9)x d9 dx 

'0 ^ Jo Jo 

P projected area of helix (at infinity) 

5 control surface (at Infinity) 
a volume of wake region 

p density of fluid 

V interference velocity 

v_ axial interference velocity 
z 

Vp radial interference velocity 

v^ tangential interference velocity 

0 velocity potential 

T thrust 

/l V + w \ 

H pitch of wake helix (2it — - — j 

6 axial energy-loss factor 
€p radial energy-loss factor 

€^ tangential energy-loss factor 

E energy loss in wake 

/ TV \ 

T] efficiency! ] 

\tv + eJ 

c specific loading factor referred to wake at 

infinity / —~\ 
\^PpV^y 

a apparent induced displacement velocity at 

propeller disk 
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MASS COEFFICIENT K 



In reference 1 the concept of a mass coefficient K 
was introduced. By definition 




whe re 

K(x) - 



PPoo 



2Tr(V w)v/ 



and X is the nondimensional radius of the wake. In 
more general terms to include also dual-rotating pro- 
pellers, K may be defined by 

Where K(x,0) is a function of both the radius x and 

the angle 9 or the time t. The coefficient k is 

thus the mean value of the circulation factor K(x,9) 

over the area of the wake cross section. 

It is shown In the follov/ing discussion that the 
momentum 

V da 



contained in the space a enclosed between two infinite 
planes perpendicular to the axis at a distance between 
them equal to the distance between successive surfaces 
of discontinuity is equal to 

KwP 

The designation "mass coefficient" originates from this . 
relation. The mass of air set in motion with a veloc- 
ity w is of the cross section kP or, if the coluimn F 
is considered to be in motion, it will attain the mean 
velocity Kw; hence the term "mass coefficient" is used 
for K. Reference 1 gives the circulation function K(x) 
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or K(x, G) and the mass coefficient k for all signifi- 
cant cases of sinrle- and dual-rotating propellers. It 
may be seen later" that the nass coefficient k times w 
is not exactly identical with the thrust coeffi- 
cient -Co because k refers to a certain momontuin and 

C 

not to the thrust* 



GEJTERAL CONSIDERATIONS Q? PROPELLER .VITH 
IDEAL CIROULATI'^N DISTRIBUTION . 



Expressions will be given for the thrust and energy 
loss for both single- and dual-rotating propelle.s. The 
discus£iion is restricted to the case of the ideal circu- 
lation distribution. In this case, the surface of 
discontinuity moves backv/ard as a rigid surface at a 
constant rate of :notion w. The equation of motion may 
be v/ritten in the form 

1 o 

h ^ o 2 ot 

whore the subscrirt o refers to the condition at 
infinity with the medium at rest. Because of the 
stipulation that the :'3ntire field moves backward as a 
rigid body, the relation 

^ = f(z-;7t,x,8) 
exists and, consequently, 

dt hz 
it follovvs that 
oGf 

— ^ = -v/v„ 

(3t ^ 

and the equation of motion for this type of rigid- 
pattern displacement flow i: , in general, 

P " Po \^"^^ ~ ^^■'■^'z 
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CALCULATION 0^ THRUST 



In calculating the thrust of the nro^eller, it is 
convenient to e^oloy an imaginary control surface, \A^hich 
encloses the nroiDeller but is infinitely distant from 
it. The control surface may be chosen as a cube v;ith 
infinitely long sides and with the propeller at the 
center and the wake directed perpendicular to one wall S, 
crossing it in the middle and extending infinitely far 
beyond or outside the control surface. Let the center 
of the wake be the z-axis.. By methods of classical 
mechanics^ the instantaneous thrust is obtained as 



= J[p - + p(V + vJv^]dS 



Introducing equation (1) transforms this integral to 



/[■ 



(V -f w)v^ + v^S . iv2 



2 



dS 



Since the thrust may vary with time as is the case 
for the ideal dual-rotating propeller, an integration must 
also be performed with respect to time. This integration 
results in the expression 



dS dt 




V2 are 



Since 0 and therefore the velocities v and 
functions of z - wt, this integral may be obtained as 
a volume integral taken over a volume of an infinite 
cross section S and a length along the axis z equal 
to the distance between successive vortex sheets. This 
distance is 



H V + w 

P OOP 



where p is the number of blades. The integral thus 
becomes 



T = ~p I [^(V + w)v^ - iv^ + v^sj 



da 



■H 

p uo 
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where the surface Integral is taken over one turn of the 
vortex surface with dS as the projection on the sur- 
face S perpendicular to the axis z. With 

r, 2Tr( V + w )w 

r ^ — K{x,9) 

PO) 

the following relation is obtained: 

~ J v^ da = wj K(x,9) dS 

This relation may be transformed by the introduction of 
the mass coefficient k (see reference 1), which is 
defined as the mean value of K(x,9) over the projected 
wake area 

K - - K(x, 9) dS 



Then 



P 

For the second integral occurring in the expression for 
the thrust, a similar treatment yields 

/ v2 da = KW^F 
±^ Jo 

Finally, by definition of a quantity €, which may be 
recognized as the axial energy-loss factor, the third 
integral is 
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This integral is obviously the expression for twice the 
axial energy loss contained in the volume c between 

successive vortex sheets. Since ^z^<^^ 9 evident 



that €<K and 



<1. 



By use of these three expressions, the thrust mo^j 
be written in the sirm^le form 



- Pp 
= pFKw 



|^(V -f w)KW + - iKW^I 



V + w 



(2) 



CALCULATION OF EITSRGY LOSS IN WAI<E 



By methods of classical mechanics, it can he shown 
that the expression for the instantaneous energy loss in 
the wake Is 



E = 



J - Po + |pv2^V2 + |pv2v 



dS 



By use of relation (1) this expression becomes 

0 



E - p I I V2 



2w + iv^^dS 



By integrating over time and transferring to a volume 
integral as was done for the thrust, the exr)ression is 
changed to 




E = P ^ r fv^^w + 
pHJa 



^v^v^da 



Replacing the integrals j da and v2 da 



before gives finally 



E = pFKw4 -^w + 



(5) 
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EFFICIENCY 



With the thrust and the energy loss known from 
equations (2) and (3), the efficiency defined as 
TV 



TV + E 



is given by 



p FKw M 



n = 



pPKw 



V +• w 



or 



V + wi 



11 = 



(I • 0 



V 



(V + w) 



With the Introduction of a nondimensional quantity w - w/v, 
the efficiency becomes finally 



1 + w - 



(1 + w) 



(li.) 



which is the exact expression for the ideal efficiency of 
the heavily loaded single- or dual-rotating propeller. 
The efficiency is a function of the velocity ratio w/v 
infinitely far behind the propeller and is dependent on 
only one other parameter c/k, which is the ratio of 
the axial loss to the total loss. 

By introducing the specific loading factor Cs in 
equation (2), the following expression is obtained: 



1 ^ T 

p ^s " p 
^ FpV^ 



= KW 



1 ^ + 0 
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Substituting the quantity c in equation (I4.) gives 



s 

^ 

(5) 



where 



This formula sbows the efficiency as a function primarily 
of the parameter Cq/k v/ith a slight dependence on € /k , 
The efficiencies from the formulas given are plotted as 
functions of w and Cg/nc in figures 1 and 2, 
respectively; numerical values of the efficiencies are 
listed in tables I and II. 

It is of some interest to give the exact expression 
for r) as a function of w and of CsA the form of 

infinite power series. These series are 

and 

The verv small dependence of on €/k, particularly 

In equation (7), Is evident. 



INDUCED \^LOCITY AT PROPELLER 

If the efficiency Is v^ritten 1.n the form 

1 
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ths quantity'- a gives the apparent induced displacement 
velocity at" the propeller disk. Frora equation (if). 



1_ c -2 
— w + — vv 

a = 



1 ^s 1 ^fsY 1.A ^ A A A/V 

^" " d V \/ 16\2 V\2 ' VV^. 



(8) 



-A c\ 

1 -H W - + - ) 

or, in power series in w, 

_ 1_ lA cV2 ^ l/l ^^^--5 ro^ 
2 2\2 2\^^ 

or, in terms of Cs/k, 

.2 .A A A/cA5 



+ ... (10) 



FINAL REMAaKS 



Exact formulas for the ideal efficiencies have heen 
developed. The efficiency is given simply as 

•0 

' -J 

It has been shovm hy series developments and graphs that 
the dependence of t) on the parameter e/K, which is 
the azial loss in terms of the total less, is very small, 
nartl;^ularly in the second formula. For this reason 
the namerlcal value of ^/^ need not be known to a 
high dogrse of accuracy. It can be shown that e/K is 
appro:."i^,ately eqval to k and, for most practical 
Durooses, this approximation is sufficient. On the other 
hand, che ^^ormulas are exact and the value of e must 
be obtained to the degree of exactness actually desired. 
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Values of € for single-rotating two-hlade propellers 
are given as an example in the appendix. 
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Appsroix 

DETExRMINATION OF AXIAL LOSS I^AGTOR € 

It is seen from the efficiency formulas that the 
axial loss ratio c/k enters as a parameter. Since the 
dependence is verj^ small, it is sufficient to know an 
anproximate value of ^-/k . It is concluded from the 
following discussion that the loss ratio c/k: is only 
slightly greater than the numerical value of the mass 
coefficient k, since this relation holds for the known 
case of an infinite number of blades and shows reasonable 
agreement also in the case of a two-blade propeller 

for ^ = ^. Until the loss ratio has been obtained by 

2 c 
direct calculation, the practice of putting ^ is 

considered satisfactorv for all purnoses. 

Axial, tangential, and radial loss factors are 
defined for single-rotating propellers by 



w^F 



2 dS 



and 



2 



r - — I 



respectively. Further, the total loss factor is given as 



de ^t ^^r 
In figure 5, — ^ — , — , and ^ are plotted 



d(x2) d(x2) d(x2; 



arainst for the case of a tv;o-blade r^ropeller 

w3th X = These TDlots were made by using the 

2 
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functions and constants given by Goldstein for the 
velocities v^, v^, and v^. The curves, upon 
integration, yield the values 

e = 0.0925 

€^ = 0.0768 

= 0.0952 

The sum of these three, 

K =: C + -f = 0,2625 

is very nearly equal to the value of k obtained 
from figure 3 in Dart I (reference 1). It is noted that 
the radial loss, which has been neglected in all previous 
discussions of this sublect, is the largest of the three 
losses . 

In the case of an infinite number of blades, the 

formulas for and — can be integrated 

d(x2) d(x2) 
explicitly to give 

c = 1 4- 



1 * \2 




^t 



o ^ log fl -f 

1 -f X" \ 



and 



c - 0 
r 



The total energy loss is then 
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The functions e^., and k are plotted against 

in fi,<^ure A plot of the function c/k is also shown 

in flsurs 1+. The value of c/^^ for a two-blade pro- 

neller with ^ which is obtained from the given 

data, is shown by a point In figure Ij.* It is noted that, 
for the case of an infinite number of blades, c/k is 
slightly greater than k; whereas, for the two-blade 

propeller with ^ = |. = 0*35 and k =^ O.2625. The 

quantity e/K may therefore be tentatively estimated 
as somewhat greater than k. 
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r] - 



(1 + w) 







w 


0 


1/100 


1/10 


1/5 


2/5 


5/5 


1 


0 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


.05 


0.9762 


0.9762 


0.9760 


0.9759 


0.9757 


o.975i+ 


0.9751 


.10 


0.951+5 


0.951+5 


0.951+1 


0.9557 


0.9528 


0.9520 


0.9501; 


.15 


o.95i+8 


0.95i+7 


0.9558 


0.9529 


0.9511 


o.929i| 


0.9265 


.20 


0.9167 


0.9165 


0.9150 


0.9155 


0.9105 


0.9077 


0.9028 



TABLE II 

PROPELLER EFFICIENCY (IN SERIES FORM) 



lAs^ 



.03x2 



i|\ K/ • 16 \K 



J 



16 







e /k 


1 

0 


1/100 


1/10 


1/5 


2/5 


5/5 


1 


0 


1 .0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


1.0000 


.1 


0.9767 


0.0767 


0.9767 


0.9767 


0.9767 


0.9766 


0.9766 


.2 


0.9565 


0.9562 


0.9561 


0.9560 


0.9559 


0.9558 


0.9557 


.5 


0.9577 


0.9576 


0.9575 


0.9570 


0.9566 


0.9562 


0.9560 


.1+ 


0.9200 


0.9199 


0.9192 


0.9186 


0.9171+ 


0.916© 


0.9160 


.5 


0.9025 


0.9022 


0.9009 


0.3995 


0.8975 


0,8958 


0.891+5 
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Figure 1.- Propeller efficiency. 
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Fig. 2 
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Fig. 3 
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Fig. 4 
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